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Abstract 

We show that a resummation model for the evolution kernel at small x creates a bridge 
between the weak and strong couplings. The resummation model embodies DGLAP and BFKL 
anomalous dimensions at leading logarithmic orders, as well as a kinematical constraint on the 
real emission part of the kernel. In the case of pure gluodynamics the strong coupling limit of 
the Pomeron intercept is consistent with the exchange of the spin-two, colorless particle. 



1 Introduction 



The high energy hmit of QCD is one of the most important aspects of strong interactions. The 
abundance of precision data from high - energy coUiders, hke HERA and Tevatron, enables one 
to confront and check the theoretical predictions in QCD. The standard framework used to make 
the predictions in QCD is that of the coUinear factorization. In this approach the hard matrix 
element is computed within the perturbative QCD. The coUinear singularities are factorized into 
the parton densities which are then evolved using the renormalization group equations. These 
equations provide a tool for computing the change of the densities of quarks and gluons at the 
variable scale. Although the parton distribution evolves, the total energy-momentum is conserved. 
This is reflected by the fact that the anomalous dimensions of the DGLAP equations vanish at 
j = u; + 1 = 2 

2Nf-/qg{uj)\^^i +799(^^)^=1 = 

and 

lgq{uj)\uj=l + lqq{uj)\u; = l = , 

order by order in the perturbation theory. Therefore parton distributions Fi{x,^'^) obey the 
momentum sum rule 

' dxxFi{x,ii^) = 1, 







where x is the fraction of the longitudinal momentum of the hadron carried by the parton, and 
H is the scale of the DGLAP equation. In the distribution Fi partons are on-shell and have zero 
transverse momentum. This by itself is already an approximation, since one does not observe free, 
on-shell quarks and in any hadronic process the target partons are always off-shell. 

An alternative approach is BFKL resummation [1] , where gluon emission diagrams are summed 
in the Regge limit s 3> The amplitude computed perturbatively in this regime has the form 

A{s,t) - , 

where ap is the Pomeron trajectory given by the integro-differential equation 

df{x,kT) _ asNc 



dlnl/x 



d^k'j,K{kT,k'rj,) f{x,k'j,) . (1) 



In ID) K{kT, k'rp) is the BFKL kernel which depends on the transverse momenta fcy. In the leading 
logarithmic approximation the Pomeron trajectory intercept equals 

ap(0) = 1 + 41n2— ^ . 

The solution to Eq. ((T]) is the unintegrated gluon distribution function /(x, fcr), which still 
contains the information about the transverse momentum of the gluon. The distribution / is for 
the off-shell gluons. 

The next-to-leading order corrections [21 [3] dominate the leading order. Consequently, several 
methods on the resummation of the perturbative series were developed, to name [l]-[9], p!0]-[15j. 
[T6] and [17]. There is, as yet unfortunately, no unique resummation procedure at small x. However, 
the major building blocks are common to all approaches, like the DGLAP anomalous dimension, 
kinematical constraints, running coupling and the momentum sum rule. 

The unique sum of the perturbation theory, if it exists, should contain the complete information 
about the anomalous dimensions. It is then legitimate to ask what is the strong coupling limit 
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of such a sum. While one does not yet have a complete solution, one can still build models and 
explore them in the strong coupling limit. In [18l [19] a resummation procedure was proposed 
and it revealed that the intercept in the limit of the infinite coupling equals one. In this paper 
we investigate the details of the behavior of the BFKL eigenvalue restricted by the kinematical 
constraint and the momentum sum rule. We show that in the limit — > oo the intercept of the 
resummed Pomeron becomes exactly one. This is consistent with the exchange of the (possibly 
massless) particle of spin two. The model proposed here explains the strong coupling limit but also 
allows to interpolate between small and large values of as ■ The importance of the exact kinematics 
for the parton distributions in Monte Carlo generators has been already emphasized [501 HI] ■ The 
systematic program of the definition of the coUinear factorization with exact kinematics is under 
way [22] . This paper solely focuses on the BFKL Pomeron. 

Our results may bear some relevance in the context of the string/gauge duality conjecture 
proposed by Maldacena [23] . Recent string theory results 2A'l (see also [25l [2^ ) showed that the 
scattering amplitudes in the Regge regime at large N and for the vanishing beta function exhibit 
a diffusion in the fifth dimension of the curved AdSs space. The kernel has the form 



with 



K(r,r',s) = ^ p-(lnr-lnr')V(4Jlns) 

V AttD In s 



io = 2-A + o(i/A) 



where A = R*/a'^, R is the radius of the AdSs space and the string tension equals l/(27rQ;'). In 
the case of the supersymmetric Yang-Mills theory A — g^N — ^nagN . Therefore the diffusion 
properties of the BFKL Pomeron are universal at all values of the strong coupling. We show that 
the resummed eigenvalue is consistent with the results of [24j . 



2 BFKL evolution and higher order corrections 

In the leading logarithmic approximation the BFKL evolution equation [T] in the momentum 
representation reads 



, (2) 



where x is the Bjorken variable, and kx, k'j, are transverse momenta of the exchanged giuons. For 
the unintegrated gluon distribution function /(cc, kr) the solution to the Eq. ^ behaves as 

/ ~ X--- , (3) 

where 

EE ap(0)-l = 41n2^^ = Xo(7 = l/2), (4) 

TT 

and the kernel eigenvalue xo has the form 

Xo(7) = -27b - -(/-(t) - ■0(1 - 7) , (5) 
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where ip is the digamma function and 7 is the Melhn conjugate to Ink^. The next-to-leading 
corrections to the BFKL equation are known to be very large [21 [3] . The resummation procedure 
yields a stable result for the intercept. The following features are common to all methods of 
resummation 

• Full DGLAP anomalous dimension (at least in the leading order of perturbation theory). 

• Kinematical constraint imposed onto the gluon emission term (explained below). 

• Running coupling constant. 

• Momentum sum rule. 

It has been shown, 12J that the Pomeron intercept computed from the renormalization group- 
improved BFKL equation with modifications mentioned above, is much smaller than the leading 
logarithmic value. The resummed value increases with the coupling constant. The correction 
which comes from the expansion of the resummed result coincides with the NLLx [51 [3] calculation 
down to a couple of percent. 

The kinematical constraint [571 [551 dS] (for Monte Carlo results within dipole picture see [501131] ) 
is imposed onto the BFKL equation as a requirement that the virtualitics of the exchanged gluons 
are dominated by the transverse components of the momenta. The form of the kinematical con- 
straint depends on the scales in the scattering process [32] . For DIS reaction which is characterized 
by a large difference of scales of the scattered particles (where virtuality of the photon is much 
larger than the scale on the proton side ~ ^qcd) kinematical constraint imposed onto the real 
emission term in ([5]) approximately equals 

k'S <'^. (6) 



In the case of the scattering of objects with similar scales (as for example in 7*7* scattering 
with photons having both large and comparable virtualities) , the kinematical constraint has a 
symmetrical form 

e[kT/z - k'rp] e[k'rr/z - kr] . (7) 

The exact form of the constraint, as it appeared in [27l[28|, actually is 

1 — z 

< fcy , qt = kr — k'rp . (8) 

z 

In the limit of small z we neglect it in the numerator of ([8]) and set k'^ ~ q^. The comparison 
between the two alternative versions of the constraint ([6]) and ([8]) is relegated to the final section 
of this paper. 

More exact analysis of the kinematics and the constraint can be best illustrated by the following 
simple example. In Fig. [1] we show a 2 — + 3 process with one particle emission. The details of the 
matrix element are irrelevant for the purpose of this discussion. The incoming momenta of the 
particles satisfy the condition — — ^ ^- The general expression for the phase space of 
such a diagram is 

j dPS,^ J ^^6m~k,f)Siip,^k,f)6{ip, + k,)'). (9) 
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Figure 1: Elementary diagram for the one particle emission in the high energy scattering. 

We decompose the exchanged particle momenta in terms of Sudakov variables 
ki = {aip^ , PiP2 J kix). The phase space can be rewritten as 



dPS. = 



4(27r) 



daidPicfkiT da2dP2d^k2T ■ 6 (s(ai - a2)(/3i - P2) - [kir - fc2T)^) 



5{-s{l - ai)f3i - klj.) 6{s{l + (32)a2 - k^ 



2Tl 



(10) 



In the limit of the multi-Regge kinematics 

0<Q!2<ai<l, 0< < I/32I < 1 
which leads to the following approximated form for the phase space 
£.2 



dPS. 



4(27r)5 



daida2d(3id(32 d^kix d^k2T 



■ 5[sail32 + [kiT - k2TY) 5{-sl3i - kfj.) 5{sa2 - . 

Furhermore, in the high energy limit the transverse momenta are approximated as k\rp ~ k2x 
which leads to 



(11) 

^ ^2 



/I /" ^ d(y. 
dPSz = / da2 5{sa2-klj,)d^kiTd^k2T 



1 



4(27r) 



d^kiT d^k2T ■ (12) 



Phase space becomes completely factorized into transverse and longitudinal components. This is 
also transparent in Eq. ([2|), where we identify x = /s, ai = xj z. From the point of view of the 
leading logarthmic accuracy, i.e. (as In s/m^)", the choice of w? is irrelevant since it has to be 
much smaller than the very large energy s. This, however, leads to a mismatch in the kinematics for 
the outgoing particle. The approximated phase space is now completely factorized, as in Eq. (jl2p . 
and the kx integration is unrestricted. This can give an unphysical result for the emitted gluon 
since it can be off-shell 



2g+g- 



1 

9t 



< if A:t 7^ fcy 



(13) 



The kinematical constraint ([5]) guarantees that in the limit z 1 the emitted gluons are soft and 
on-shell ( ~ '^Q^^T — (Zt — 0)- Note that the condition ([5]) constrains only the phase space. 
The matrix element is kept in the same, high-energy approximation. 
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3 BFKL with kinematical constraint 



Let us now explore the details of the effects of the symmetrical kinematical constraint ([7]) imposed 
on the BFKL kernel. Throughout this analysis we will keep the strong coupling fixed. The modified 
equation read^ 

/(x,4) = /W(a;,4) + 

dz f k^dk'^ \ /(f , fc^2)e(/c^ _ kTz)Q(kTlz - k'j) - /(f , /cl) 



t'2 1 |l.'2 _ 1,2 I 



[4fc^ + 4]i/2 



(14) 



where we introduced standard notation as = a^iVc/Tr. We perform the following change of the 
variables 

77 = ftslni , 

z 

^ = Ssln- , (15) 

X 

after which the equation ((M]) can be rewritten as 

+ f^dr^ f ^ ^^'^ ' ^' ^^'^ - ^exp(-r;/«,)] e[A:T exp(77/a,) - fc^] - f{ii - fcf,) 



y2 1 I 1./2 _ iL2 



/(>?-e,fc|.) 

[4fc^ + 4]i/2 



(16) 



The LO BFKL is recovered when the theta functions are set to one and its solution can be written 
in a form of series in ^ 

i 

which is the leading logarithmic series in aslnl/x. Clearly, the non-perturbative nature of the 
kinematical constraint is evident because of the appearance of the exponential factors exp(±ry/Q;s). 
It is non-perturbative in a sense that the leading logarithmic expansion in rj is accompanied by 
the arbitrary powers in Ois- The resummed Eq. [16] contains all powers of ols as well as ^ and the 
solution is a fimction of three variables f{£,,as,k^), or f(x,as,k^). 

For the fixed values of 77 variable the O functions have the following behavior 

e[A:^-A:Texp(— ^)]e[fcrexp(^)-fc^] 1, (17) 

as as 

e[k'j,-kTeM- — )]<^[keM — )-k'T] "'-^ —5(kT~k'rr). (18) 

as as as 

The constraint breaks down the complete factorization between the longitudinal and transverse 
components of the momenta in the evolution. The evolution equation cannot be written any longer 
in a simple form of the differential equation with the factorized In 1 /x dependence as in Eq. [TJ 



^We write here f{x,ki,) instead of f{x,kj-) to emphasize that the angular dependence has been integrated out. 
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The Mellin transform of the kernel in Eq. (|14p reads 



1 dz 



- 1 



I krp krp I 



[4fc^4_^fc|.]l/2 



1 dz 



j/2 



dii u''' — 1 
u 1 — li 



duu'^ — I 



du 1 



-27B-i^(7 + |)-VXl-7+|) • (19) 




Figure 2: The behaviour of the eigenvalue function xil:^) i Eq. [121 at 7 = 1/2 and around u = \. 

The effect of the kinematical constraint is to shift the arguments of the digamma functions. The 
arguments 7+u;/2 and 1—7+^/2 can now be equal or larger than 1 in the domain < 7 < 1, w > 0. 
Therefore, both BFKL branches, '0(1) ~ "0(7 + w/2) and ■0(1) — ^"(1 ~ 7 + w/2), have zeros in this 
regime. This results in a solution x (70, ciJo(7o)) — with < 70 < 1. In particular, the eigenvalue 
x(7,w) (I19p vanishes when 7^1/2 and lo ^ 1 

x(7 = l/2,co-l) = 0, (20) 

which is illustrated in Fig. [2] The solution to the eigenvalue equation 

u} = OLsxil,^) , (21) 

can be symbolically written as [13j 

u; = x°''(7,a.) • (22) 

It is difhcult to obtain the solution analytically, however we can easily solve numerically for 
In Fig. [3] we illustrate this solution as a function of 7 for different values of a^. For increasing 
values of as, the minimum of x''^ at 7 = 1/2 approaches the limit of 1 as is expected from Eq. 1201 
The shape of the function changes very little with the increase of cxs- The derivative of xij,'-^) at 
u) = 1 equals 

dx(7-l/2,^) i TT^ 

dl^ ^-=' = "T ■ 
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Therefore, we can write an expansion of the eigenvalue (j2ip around cu = I 



X(7 = l/2,c^) ^ 



(24) 



however, the shifted kernel eigenvalue is not sufficient to reproduce the full NLLx BFKL result. 
Motivated by the resummation procedure proposed in |12[ |13| let us consider the resummed model 
for the eigenvalue equation in the form 



1 



-+A{lj) 

id 



The expression in the square brackets is the LO DGLAP anomalous dimension 



1 



UJ 



(25) 



(26) 



Xeff (Y, Os ) 

2.5 




-0.5 



0.5 



1 . 5 



Figure 3: The solution to the eigenvalue equation (|22p in the case of the LO BFKL, with kinematical 
constraint in the form of the function x°*(7, as), for four different values of the strong coupling: 
cxs = 0.2, 1.0, 5.0, 20.0, from the bottom to the top curve. 



A{oj) = - 



1 



1 



1 oj + 2 
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with 



and 



(27) 



(28) 



7 



(0) 



= 



Ail) = 1 . (29) 

The last condition, expressed by formulae (|29p . is the momentum sum rule in the absence of 
quarks Nf = 0. Therefore, the resummed model (I^S)) contains DGLAP anomalous dimension at 
to leading order in ofglnQ^, BFKL eigenvalue at the leading order in Inl/x and the kinematical 
constraint. We again stress, that the coupling constant is fixed in this analysis. In particular, this 
enables us to use simple expressions for the eigenvalue conditions like Eq. (PTjl . We expand the 
DGLAP anomalous dimension 7^''^ (w) around lj = 1 and get 



7 



(0) 



^ '^6 144 



= 1.1935 (1 -tj) 



(30) 
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Using (gH) and wc can thus write (HH) 



' = (l--)'^(^-S) - l-9633(l-c.)^ (31) 



Us ' 6 V 6 144 

Solving for lu we obtain 

CJ = 1 r= 

where 



(32) 



CO ^ l/yy^y-Yiij - 0-509346. (33) 

The result ([5^ should be compared with the calculation of the first correction to the intercept 
of the graviton presented in [25 , 

j= 1+^^ = 2 Co = I/tt, 



and the SUSY result [Iff,. The coefficient co obtained from the model ((25)) is different, due to 
different details of the resummation model and the fact that here we are considering fixed coupling 
QCD, whereas Refs. [Tn| apply to the SUSY YM case. However, the overall agreement is 
satisfactory. 

Xeff (Y,as) 
2.5r 



2 



1.5- 




-0.5 0.5 1 1.5 

Figure 4: The solution to the eigenvalue equation in the case of the LO BFKL with kinematical 
constraint and the DGLAP anomalous dimension in the form of the function x°^(7, ag) for the 
five different values of the strong coupling as = 0.2, 1.0, 5.0, 20.0, 100.0 from the bottom to the 
top curve. 

The model (l?5|) provides not only the correct limits at small and large values of the coupling 
but also smoothly interpolates between these limiting regions. From (j3ip it is clear that the l/-\/a7 
behaviour is a results of the double zero of the model ([25]) at u; ~ 1. The numerical solution to the 
equation ([25]) in the form of H 

^ = X°*^(7,as) , 
■^We put an additional tilde'to distinguish it from the result of Eg. 1221 
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is plotted in Fig[3]where we show it for various values of the coupling constant. The fixed points in 
the (7, w) plane, namely (—1/2, 1) and (3/2, 1), result from the superposition of the zero lo = 1 and 

the simple poles at 7 + a;/2 and 1 — 7 + cj/2. The second derivative ^ d^'^'"'^ It^i/z appreciably 
changes in this case, exhibiting l/ag behaviour at large as- This is different from the AdS/CFT 
correspondence result which as well produces for the second derivative. We think that this 

is an artifact of the simplified model which we are considering. The overall qualitative behavior 
of is nonetheless consistent with the expectation from string theory side, in particular when 
we compare with the Fig. 2 in Ref. f24;. The effective eigenvalue has two fixed points (the points 
which do not change when the coupling is varied) and it can be approximated by the parabola 
with a minimum at 7 = 1/2. When the coupling is increased, the second derivative decreases. This 
continues until the limit of ag = 00 is reached, and the effective eigenvalue tends to a constant 
equal to one. 

When we compare the solutions presented in Fig. [3] and Fig.lH we see a qualitatively different 
behavior of the effective eigenvalue function. In Fig. [3] the function x'^^ is 1 only at a single point 
in 7 = 1/2. In Fig. [4] the effective eigenvalue bmction has a constant limit when — > cx) for 
all values of 7, which is the consequence of the fact that 1 + ujA{uj) ^ {1 — uj) at w = 1. 

From x°^(7j'5s) one can compute the dependence of the Pomeron intercept uip — x'^*''(l/2, a^) , 
on the coupling which we illustrate in Fig. O compared with the first order expansions around 
cis ~ (leading logarithmic BFKL [T|) and ~ 00 (the string theory in the curved background 
[24]). 




agN/TT 



Figure 5: The strong coupling dependence of the Pomeron intercept in the case of fixed coupling. 
Comparison of the predictions of the resummed model with LO DGLAP (|25|) (solid line) to the first 
order results in QCD [1] (dashed line at small couplings ) and strings in curved background (dotted 
line at large coupling) [24^ . Note the logarithmic horizontal axis for the coupling as = "'^^ ■ 

It is also interesting to examine the convergence properties of the resummed model at small 
and large values of cts- The coefficients of the expansion for the intercept cup — x°^{l/2, as) are 
collected in Table 1. The left column contains the first five coefficients of the series for the Pomeron 
intercept at small values of as 

wp = ^ a„ a" , 

n>l 
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and the right cohimn contains the first five coefficients of the series at large values of as 

ri>l 

The coefficients in the small coupling expansion grow with increasing order as a„+i/a„ ~ — 10 
signalling the possible lack of convergence of this model at 0. At large values of the coupling 

the coefficients behave at least as 6„+i/5„ ^ 1/n. 



an 


bn 


2.77 


0.713685 


-20.72 


0.300767 


226.69 


0.085747 


-3000.73 


0.015120 


44267.30 


0.000863 



Table 1: Numerical values of the expansion coefficients of the intercept at small (left column) and 
large (right column) values of the strong coupling constant as- 



4 Improved kinematical constraint 



It is interesting to investigate the changes to the evolution kernel when imposing more stringent 
kinematical constraint. In [28j it was shown that the more accurate version of the kinematical 
constraint has actually the form 



1 - z 



(34) 



It reduces to Eq. ^ when z <C 1. Let us first consider the coUinear approximation to the BFKL 
kernel 



/(a.,4) = /(°)(x,/c2) + 





' /■'^^ dk'^ 


F krpdikrp 


Ix z 


J k!^ 





f{x/z,k'^) 



(35) 



and keep the (1 — z) term in the numerator of the kinematical constraint replacing however ~ fc^. 
This constraint also divides the domains of the integration over the longitudinal variable z. There 
are three distinct regions. For large values ofz> 1/2, (1 — z)/z< 1 and k!^ < k^ we have 

= -A/2(l + 7,^-7) = ^-^2^^i(l + 7,l + ^;2 + 7;-l) , (36) 
7 7(1 + 7) 



and for the region where k!^ > k^ and z < 1/2, (1 — z)/z > 1 we obtain 

dz 

lo z 



coll.u f \ 

X.<i/2(7,^) 



/ = If ) 6(4(1 -.)/.- fc^^) 6(1/2 -z) 



T \ "-T / 

- A/2(l-7 + ^,7- 1) = ^ , . 2j^i(l-7 + ^,-^;2-7 + ^;-l) , (37) 
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and finally z < 1/2, fc^ < kr 
The eigenvalue equation for ([?5]l is then 

1 - I coll.l ( \ , coll.u f \ . colli / \\ 

1 = (x,>i/2(7, ^) + X.<i/2(7, ^) + X^<i/2(7, j ■ 

The effect of the improved kinematical constraint on the coUinear model results in the shift of 
the pole at7 = l— !-7 = l+ cj. Additionaly the regularized beta functions generate the additional 
poles at 7 = — 1, —2, —3, . . . from Xz>ili(n^ ^^i^ 7 = 2 + oj, 3 + tj, . . . from Xz<i/2(77 These 
singularities are located in the same positions as those in the BFKL eigenvalue ([TO]) with an 
asymmetric shift of the 7 poles. This has to be compared with the following expression 

lx^°"(7,..) ^ + ) , (38) 

w ^V7 1 — 7 + ti-'/ 

in which the small-z version of the constraint fc'^ < j z imposed onto the coUinear model leads 
to the shift of one of the existing poles in 7 space but does not introduce any new singularities. 

We next impose the constraint < k^{l — z)/ z in the BFKL case 

f{x,kT) = .f^"\x,kT) + 

In this case we perform the diagonalization of the kernel numerically. The eigenvalue xili ^) a 
function of 7 for different values of uj is shown in Fig. [51 where it is compared with the expression 
(fT^ for the case of the LLx BFKL with the approximated kinematical constraint. The eigenvalue 
is further lowered with respect to but the overall differences are moderate. The most affected 
region is the one in the vicinity of 7 ^ 0. This has to be understood as an effect of the large value 
of z in the cutoff (1 — 2;)/z < 1 for z > 1/2, and therefore the cutoff becomes more stringent than 
the previous one in the coUinear region where k'^ < k'^. 



2 



qi 



(39) 



5 Conclusions 

The BFKL equation in the leading logarithmic approximation suffers from the problem of the 
violation of the energy - momentum conservation. The kinematical constraint reduces the phase 
space allowed for the the real emissions and is responsible for the resummation of the subleading 
corrections. The kernel eigenvalue in that case is zero when oj — 1 and 7 — 1/2. This means that 
in the limit cis 00 the intercept becomes equal to 1, which can be interpreted as an exchange of 
the color- singlet object with spin 2. When the momentum sum rule is additionally imposed, the 
effective eigenvalue becomes constant and equal to one in the limit of the infinite strong coupling. 
Such model then can provide an interpolation between the small and large values of the strong 
coupling constant. The first correction to the intercept of the graviton, of the form l/^/al can 
be also computed from this model. The more accurate form of the kinematical constraint further 
reduces the eigenvalue. However, the overall qualitative behavior remains unchagcd and the model 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

y 



Figure 6: The eigenvalue of the BFKL kernel with the kinematical constraint of the form ([6]) 
dashed lines, and ([8|) solid lines. Calculation performed for three different values of the Mellin 
variable lo = 0.1, 0.5, 1.0. 

still has to be supplemented by the requirement that at w = 1 the eigenvalue vanishes for all values 
of the Mellin variable 7. 

This result is valid in the case of one Pomeron exchange. We stress that the result for the 
intercept will be renormalized when the multi-Pomeron exchanges or the saturation corrections 
are also incorporated. These corrections will certainly become important in the limit of the very 
large coupling. Also, the running of the coupling and the presence of quarks in the evolution can 
further change the results in a quantitative way. We leave these questions for further studies. 
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